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Abstract 

We study the low-energy effective theory in N = 2 SU(N C ) supersymmetric QCD 
with Nf < 2N C fundamental hypermultiplets in the Coulomb branch by microscopic 
and exact approaches. We calculate the one-instanton correction to the modulus u = 
(ITtA 2 ) from microscopic instanton calculation. We also study the one-instanton 
corrections from the exact solutions for N c = 3 with massless hypermultiplets. They 
agree with each other except for Nf = 2N C — 2 and 2N C cases. These differences 
come from possible ambiguities in the constructions of the exact solutions or the 
definitions of the operators in the microscopic theories. 
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The prepotential of the low energy effective theories of the N = 2 supersymmetric 
gauge theories in the Coulomb phase have been obtained exactly by using the holomorphy 
and the duality arguments The prepotential receives a non-perturbative instanton 
correction in the semi-classical region [gj. On the other hand, the microscopic instanton 
calculations in supersymmetric gauge theories || f|, || give reliable results when the semi- 
classical approximation is valid. Thus the comparison between these approaches provides 
a non-trivial check on the method of the microscopic instanton calculus as well as the 
assumptions in the derivation of the exact solutions. 

Such comparisons have been studied in the case of N = 2 supersymmetric Yang- 
Mills theories and 577(2) supersymmetric QCDs (SQCDs) |6|-|ll0fl. The results have been 
consistent with the exact solutions so far, while some discrepancies have been reported 



in the SU(2) SQCDs with Nf = 3,4 [10]. These discrepancies do not seem to contradict 
with the assumptions in the derivations of the exact solutions in the sense that they could 
come from the possible ambiguities in the exact solutions or the definitions of the quantum 
observables in the microscopic theories. 

In this letter we study the one-instanton correction to the prepotential in the N = 2 
SU(N C ) SQCD with Nf < 2N C fundamental hypermultiplets both from the microscopic 
and the exact viewpoints. In the case of N = 2 SU(2) SQCD, the contributions from the 
odd numbers of instantons vanish due to the anomalous Z 2 symmetry since the funda- 
mental representation of SU(2) is psudoreal But, for N c > 3, one can expect that all 
the instanton corrections appear in general. 

We introduce an N = 1 chiral multiplet = (A, ip) in the adjoint representation and 
an N = 1 vector multiplet W a = {v fil A), which form an N = 2 vector multiplet. The 
N = 1 chiral multiplets Qi = Vw) and Qi = (qi,ip m i) {i = l,---Nf) form the Nf 
N = 2 matter hypermultiplets in the fundamental representation. The microscopic N = 2 
Lagrangian for N = 2 SQCD is given by 

£ = 2 J d A 9Ti (tfe- 29V <t>e 29V ) + ^ (/ d 2 6W a W a + h.c.) 

N f / N f 

+ Jd 4 6j2 {Qle~ 29V Qk + Q k e 29V Q{) + [iV2gj d 2 6 ]T QM h + h.c. 
fc=i V fc=i 
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+ ( / XI m kQkQk + h.c. 



, :i; 

k=l 

where g is the gauge coupling constant and the trace is taken in the fundamental represen- 
tation. Here the color indices are suppressed. We will examine the euclidean lagrangian 
of ([!]) in terms of the component fields in the Wess-Zumino gauge. 

The Coulomb branch of this theory is parameterized by the expectation values of the 
adjoint scalar vacuum expectation values Aq^ = a^l and A\l = a,i5{. For generic values 
of a« and a; the non-abelian gauge symmetry completely breaks down to the U(l) Nc ~ 1 , 
and the system is in the Coulomb phase. 

Since the holomorphy argument [|[ [□]] for the gauge coupling g shows that the calcula- 
tion in the region g -C 1 is enough to obtain reliable results, we will perform microscopic 
instanton calculation in the lowest order of g. In this approximation, the equation of 
motion of the gauge field D^G^ = has the instanton solutions [12 1. Their bosonic zero 
modes of the one-instanton solutions are the instanton location xq in the euclidean space, 
the instanton size p and the location in the color space. The integration over the location 
in the color space is defined by the integration over the minimal embedding of the sub- 
group SU (2), where the one-instanton configuration resides, into the gauge group SU (N c ) 
||13|| . The generators of the minimally embedded SU(2) subgroup can be characterized by 
Q j; J a Q, where Q G SU(N C ), and J a are the generators of the SU(2) subgroup obtained 
by the upper-left-hand corner embedding of the two-dimensional representation of SU(2) 
into the A^ c -dimensional representation of SU(N C ) Hence the integration in the color 
space is performed by sweeping Q in V(N C ) = SU(N C )/SU(2) x U(l) x SU(N c -2), where 
the U(l) x SU(N C — 2) is the stability group of the embedding and the additional SU{2) 
|14| comes from the fact that we are interested only in the observables symmetric under 
the space rotation. By the global gauge transformation, the group integration can be per- 
formed by rotating the scalar vacuum expectation values (A) = QAqQ^ , (A^) = QAqQ" 1 , 
while the instanton configuration is fixed at the upper-left-hand corner 0. 

The equations of motion of the adjoint fermions are given by r^D^ip = r~D^\ = 
in the lowest order of g. From the index theorem, each of these equations has 2N C zero 
modes. Similarly, each of the lowest order equations of motions of the matter fermions 
T~D^p mi = T~D^ip mi = has one zero mode. We use £, (, r\ and f] to label the zero-modes 



of A, if), ip m and if) m , respectively. Then, under an appropriate choice of the normalization 
of these fermionic zero modes[|, the integration measure of the N = 2 one-instanton zero 
modes is given by |T3|, [L5| 

210^2^+2^ -4^ f d 4 XQ f°° d 4iV c -5 f dQ f ^jUt^N, d N % 
J JO JV(Nc) J 

A^/^exp(-^, (2) 

where /i is the Pauli-Villars regulator and b± = 2N C — Nf is the one-loop coefficient of the 
beta function. 

Out of the 4N C zero-modes of the adjoint fermions, four are the supersymmetric zero- 
modes £g S , (s S obtained by the supersymmetry transformations of the one-instanton 
configuration of the gauge field. Another four are the superconformal zero-modes 
(g C obtained by the superconformal transformations. The fermionic zero-modes other 
than the supersymmetric ones, say «f = (£ sc , f a , f a ) and C' = (Csc, L, Co) (a = 3, ■ • • , N c ), 
cease to be zero-modes by taking into account the mass terms and the Yukawa terms || || . 
One of the lowest order contributions may be obtained by substituting the solution of the 
scalar field equation of motion in the lowest order D 2 A\x) = with the asymptotic value 
(A^), and the fermionic zero-modes into the Yukawa coupling term g J d 4 xTr (ijj[A^ , A]). 
Thus we obtain a bilinear term £'gM((A^})C [0 with 

gM((A^)) = ig (A)® ~^f^ c -2 + V2(A^ 

[ {£{A t } my _i^ /iVc _ 2 + v ^ ((A t ) (4) )i o 

(3) 

where e and In c -2 are a two by two antisymmetric tensor with e 12 = 1 and an N c — 2 
by N c — 2 identity matrix, respectively. Here we have divided the row and column of the 
scalar field into the following blocks; 

A ( Ail) A(2) \ 

where A&\ A^ and are 2 x 2, 2 x (N c - 2), (JV C - 2) x 2 and (N c - 2) x (JV C - 2) 
matrices, respectively, and is the traceless part of (A^)^. Another contribution of 

1 The normalizations of the zero-modes are taken to be unity under the norm J d 4 x2Tr((j)^(j)) and 
J d 4 x4>'(j) for adjoint and fundamental fields, respectively. 
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order g is the term 



N 



f 



E^^((^) (1) )+^j» (5) 



which comes from the Yukawa term gip m Aif) m . Here we have added also the contributions 
from the mass terms of the matter fermionsQ. 

The mass terms among the zero- modes (|3]), ([]) has a weight ^fg per one fermionic 
zero-mode. Therefore there exist other contributions with the same order in the case of 
SQCD. In fact, one needs to introduce the terms with four fermionic zero-modes of order 
g 2 . One of such contributions comes from the Yukawa terms gipmXq* and gqipipm mediated 
by the propagator of q. Using the scalar propagators in the instanton backgrounds |16| , 
we obtain 

12ttV ^ ^ 3 32ttV ^ 

The contribution A q S mediated by the propagator of q can be obtained in a similar 
manner. Due to SU{2) R symmetry, the sum A q S + A q S is simply given by replacing £ a £ a 
by £aCa + £aCa m the first term in (|6]), while the second term vanishes. 

The other contribution of order g 2 comes from the Yukawa terms <jTr([A, ij)]A ') and 
gip m Aip m mediated by the propagator of A, and the result is 

g 2 N f Nc 

= ~ oa 2 2 X VkVk J^iZaCa + £aCa)- (7) 
Z47F P k=l a=3 

Hence, including 87r 2 p 2 / from the contribution of the kinetic term of A, we obtain the 
classical action of the instanton configuration as 

N f i \ 2 N f N c 

S = 8n 2 p 2 f - g?MC - £ ^=Tr((A)V) + m k fj kVk - -/^ £ fj kVk J^iUa + Ua), (8) 

fc=l \V-6 / A n P k=l a=3 

where 

f((A), (At)) = Tr (<At)W<A>« + ^(A)( 2 )(At)( 2 ) + l(A^(Af^ . (9) 

The supersymmetric zero-modes must be canceled by the insertion of appropriate 
operators. An approach to do this is to consider the four fermion correlation function 



2 We treat the mass terms perturbatively. See || for another treatment. 
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of the classically massless modes of ip an d A [Q]. But it turns out that this approach is 
not so convenient for the cases 2N C — 2 < Nf < 2N C , because it detects only a certain 
combination of the fourth derivatives of the prepotential of the effective theory and this 
is always zero for the massless 2N C — 2 < Nf < 2N C cases from the dimensional analysis. 
Thus we insert the modulus u = (^TrA 2 ) || [L0|, which has the following direct relation 



to the prepotential of the effective field theory for the massless cases |17 |; 

ibiv 3T Nc BT 

^ = - A M = S«.g-^- do) 

1=1 1 

Here the prepotential has the following expansion in the weak coupling region for Nf < 
2N C : 

Ha) = T j E f E(«i - ^ lo S - IT ^ ^ bg tf) ~h ^ 

i=l \*<? i=l / n=l 

(11) 

where the first and the second terms are the classical and the one-loop parts, respectively, 
and the last ones are the instanton corrections. The logarithmic parts contribute to the 
classical part of u. Further corrections to u come purely from the instanton effects: u = 
\ Y^=i a? +ES£=i u k A blk , where u k = kT k . On the other hand, one can obtain the modulus 
u — \ Ya=\ a i + ^2h=i u T st '-^ b d k from the instanton calculation. In the massive case, u k nsL 
depends on aj and m k . In the following, we will determine the one- instanton contributions 
u\ and u % i st ' from the exact solutions and the microscopic calculation, respectively. 

The contribution of the supersymmetric zero-modes to the field A is obtained by 
solving the classical equation of the motion D 2 Ags + \f2gi[\ss, *Pss] = 0. The solution is 
given by A ss = ^ss^ssa, so we obtain f| 

9 2 ( l ?2 \ n,2 



^=-2WW-J (12) 

for the part with the supersymmetric zero-modes. 

For the massless case, after the integration over the bosonic and fermionic zero-modes, 
we obtain 

^N c ,N f < st iN c , N f ) = l N f2- b ^+ 1 A b d ) NctNf Ur t -(N c , N f ), (13) 
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where we have rescaled the field — > gcj) to match with the standard convention used in 
the exact solutions, and 



min[JV>,2JV c -4] 

U{ nsL (N c , N f ) = 2 - 5N ^ +9 - N f7i- 2Nc+4 J vn dQ Yl Nf C k T{2N c -2-k) 



x 



(Tr((A)^)) N f- k (-V2t) k det k (M((A^})) 

/((A),(At))^ c -2- fe 

Nr. 



det fe (M((At))) ee J 2 (' [YliUa + Ua) J exp (e'M((At))C) . (14) 

Firstly we will enumerate U{ nst '(N c , Nf) by estimating the structures of the poles jTJ. 
Although the integrand depends both on and a,, the holomorphy argument tells that 
U{ nst '(N c , Nf) should be a function only of the holomorphic variables a$. A pole may exist 
when the denominator of the integrand has some zeros in the integration region. This 
condition turns out that two of the coincides, because in this case (A)^ = (A)^ = 
(A)^ = is realized. Let us study the structures of the poles with the highest order. 
This comes from the k = term in the sum (|TJ]) ^. To estimate the structure of the pole 
at <2i = a,2, let us introduce an infinitesimally small parameter e by a\ — a.2 = e. Since 
f = 0(e) at n = 1, we restrict the integration region to the infinitesimally small region 
= exp(— i\faS) G V(N C ) to keep / to be 0(e). Then the nonlinear integration region 
of Q is linearized, and one obtains easily the following behavior of U[ nst '(N c , Nf): 

^^- ^JZt&^r - (15) 

The full expression should have the similar poles at a« = Oj (i ^ j), and so we obtain the 
following result up to possible gauge invariant regular terms: 

Tjinst. {N N \ EK,-(fli + aj) Nf nfcZjWi^ ~ U bi/ a i - Q fe)( Q J ~ a fc) (16) 

The gauge invariant regular term may exist only for the cases Nf = 2N c — 2, 2N C from the 

dimensional analysis. For Nf = 2N c — 2, it is a constant term, while it is a term of the form 

const. X)i=i af for the case iV/ = 2N C . These terms can not be fixed by the present method 

3 The explicit integration discussed below shows that the k — term results in the expected holomor- 
phic terms as well as some unwanted terms such as terms with logarithmic poles and non-holomorphic 
terms. These terms cancel exactly with the terms from k > 0. 
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of estimating the structures of the poles, but the explicit integrations for the cases N c = 3 
with Nf < 2N C show that the regular terms in fact vanishfj. The explicit integration is 
rather cumbersome even for the case N c = 3. This is simplified enormously by putting the 
antiholomorphic variables di to special values so that the integrand takes simple forms 
@, because the function U[ nst '(N c , Nf) should be independent of the antiholomorphic 
variables a*. Putting a\ = a,2 = 1,0.3 = — 2 and taking into account the delta functional 
contribution at Q = 1 0, we obtain (|TH) for the cases Nf < 2N C with N c = 3. 

We may also consider the massive cases. We expand with the mass terms in the 
instanton action (||). Since these terms cancel part of the matter fermionic zero- modes in 
the instanton measure, we obtain 

N f 

u\ nsL (N c ,Nf,m) = Y,t k (m)u\ nst iN c ,Nf-k), 

tk(m) = m h ' ' ' m ik- ( 17 ) 

ii<— <i& 

Now we will check the consistency of the above microscopic results fll3|) , (|l6l) and ( 17]) 
with the physical matching condition of the dynamical scales. The physical matching 
condition of the dynamical scales in the Pauli-Villars regularization scheme is given by |J 

£^Af = A^, (18) 

where denotes the dynamical scale of the original system, and A' d denotes that of the 
induced system after integrating out the heavy modes of the vector multiplets with masses 
mw,iS and the matter multiples with masses mg/s. First consider the case that some of 
the masses of the matters, say rriQ^ (i = 1, • • • , k), are very large compared to the others. 
One can show easily that, taking the limit A^ — > of ( |I?D with fixing nf=i m Q,i^- b d = ^ct 1 
will give a similar expression of fllTD with the substitution Nf — > Nf — k. Another check 
is given by the Higgs breaking (N c ,Nf) — > (A^ c — 1,0) by taking the limit b — > 00 in 
cii = a't — b(i = 1, • • • , N c — 1), cln c = (N c — 1)6. The heavy masses due to the large 
vacuum expectation value b are given by rri\y = \f2N c b and mg = \[2%b. Hence the 
physical matching condition of the scales is given by 

^d,N c -l,0 — 1 z ly c ^d,N c ,N f ■ l iy J 



'Eq. (|l6|) is correct also for N c = 2. 
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This is consistent with ([13]) and (f§ because of U X (N C1 N f ) ~ b N f- 2 N c - 2 Ui{N c - 1, 0) in 
the b — > oo limit. 

The exact solutions are determined by the hyperelliptic curve and the meromorphic 
one-form on it 0. There are some proposals with non-perturbative differences consistent 
with the symmetries of the system fl8| , [19| , [2(J. Firstly, we shall use the hyperelliptic 
18| . The hyperelliptic curve and the meromorphic one-form A for the SU(N C ) 



curves m 



QCD with Nf(< 2N C ) flavors are given as follows: 



F(xf -G(x), 

N c | 



for N f < N, 



F(x) e *) + | 2 -2 Abl ^Nf-N. x N f -N c - H .^ for Nf > N< 



G(x) = A 6i n(^ + 



rrii 



A 



i=l 

xdx ( FG' 



2my V 2G 



(20) 



The curve for the case Nf = 2N C is given by the following substitution in the above 
definitions: 

N a -2 N c 



F{x) 

Sk 

G(x) 
A 



x N ° + l(q) J] s Nc . l x L + 2- 2 L{q)Y,x l t Nc 

i=0 i=0 

(~ 1) ^ ■ • ■ ei fe , 

«i<---<ifc 
2V C 

L(g)n(^ + K?K). 

i=l 

1 xcte / FG" 



m 



F . 



l(q) 2niy \ 2G 

Here q = exp(27rrr) = exp(— 87r 2 / f?ex) an d the L(g) and l(q) are defined by 



(21) 



m 

9[mi 1112} 



40 



ho 



1(9) 







0[O 0] 4 ' v ^ 0[O 0] 

1 



exp < 27r? 

-jg^JVc-l 



(n + mi) r(n + mi) + (n + mi) r«2 



:22) 



where = r(5ij + 1) = 1, • • • , N c — 1), and and | denote the zero vector and a 
vector with one of its entries being | and the others are zeros, respectively. 



The vacuum expectation values of the scalar field dj can be written as periods of the 
one form A on the curve [ffl: 



(Ai,o) = f X, (23) 

JAi 

where Aj are the fundamental weights and A^ are the appropriate homology cycles on the 
curve. The equation (^) gives (Aj,a)'s as functions of e^'s. Inverting them, one obtains 
the modulus u = \ Y^=\ e j m terms of e^'s. For N c = 3, the explicit form of the contour 
integral (||) may be obtained by solving the Picard-Fuchs equations [[H]] with respect 
to u and v = e\e 2 e^. In the semi-classical region, the power series type solutions for Nf 
flavors (1 < Nf < 5) takes the form 

(Ai,a) = w(a Nf , f3 Nf ;x 1 ,x 2 ) + ^w('y Nf ,5 Nf ;x 1 ,x 2 ), 

(A 2 ,a) = w(a Nf ,p Nf ;x 1 ,x 2 ) - -w(j Nf ,S Nf ;xx,x 2 ), (24) 

where 

x = v l X2 = K s-N fu N f ,2^H-i) N n/^-(-v N n/2 ) (25 ) 

and a Nf = - ^znJ) » #at/ = <W/ = _ 6^v 7iv> = 4(6-^) • w(a,p;x 1 ,X2) denotes 
a power series of the form ^ mn >§d m ^ n x™ Jra x r ? 2 r ^ with rf ,o — 1- The coefficients d mn are 
determined recursively. For Nf = 6 the explicit evaluation of the contour integral would 
be effective instead of using the Picard-Fuchs equations. For massless cases we may derive 
the one-instanton correction to u for N c = 3 by the explicit enumeration of the integral, 
which is expanded in powers of e = A bl ^ 2 . The calculation goes as follows. Let e[ and 



e" denote the two branch points of the curves (20), (21) approaching to in the limit 
e — > 0. When the homology cycle in the right hand side of ( p3| ) is taken around the e[ 
and e", one obtains directly a,. The e'j and e" can be expanded in the integral powers of 
e, and we take the terms up to order e 4 . We expand the one-form A up to order e 3 , after 
the change of variable x = + ez. Performing the contour integral explicitly, we obtain 
aj in terms of e/s. Inverting the results, we obtain u — | Ylf=i e f m terms of aj's. Taking 
the term with e 2 , we obtain the one-instanton correction as 

A b ^ NfUl = 2~ 1 A b l iNj U 1 (N c ,N f ), 

TT (AT AT ^ l^i=i a i /A [ax, . . . ,di, . . . ,a Nc ) ^ 

Ui{N c , N f ) = — — — + A Nc d N 2N c -2 + B Nc d N 2Nc 2^ ( 

ii\ c {ai, . . . , a Nc ) i=1 
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A m ( ai ,...,a m ) = HK-a,) 2 ( 26 ) 

k<l 

for the cases N f < 2N C with N c = 3, where A 3 = 0. For the case N f = 2N C with iV c = 3, 
the scale parameter N should be replaced by the L(q) ~ 64g, and we obtain B 3 = — | 
by using /(g) ~ 1 — 40g. 

We have also calculated the modulus u in other curves. For the curve presented in 
|T9| , we obtain A b ^ c2Nc ~ — QAq, A 3 = — | and i? 3 = — |. For the curve [p^ , we get 
A^, 2JVc ~ -108g, A 3 = -| and 5 3 = 0. 

In order to obtain the relation between the dynamical scales and the scale parameters 
in the curves, let us consider the substitution = e^ — b (i — 1, • • • , N c — 1), e^ c — (N c —l)b 



and take the b — > oo limit in the curves (20), (21). This corresponds to the Higgs breaking 



we considered in the check of the microscopic instanton calculation, and in fact the curve 
with (N c , Nf) reduces to that with (N c — 1,0) after the rescaling and shift of x and y and 
the substitution 

7 % e -l,0 — iV c N c ,Nf ■ 

Comparing (p7|) with the physical matching condition ([19]) and using the known relation 
Ad,2,o = A.2,o @, 01 , we obtain the relation between the dynamical scales and the scale 
parameters of the curves as^ 

\h _ -Nfn-bi/2+2 \bi /r,o\ 

ly N C! N c ~ Z Z ly d,N c ,N f - 

Now let us discuss the differences between the microscopic instanton calculation and 
the exact solution. Comparing flTBD and (BBI) with using (EST) , the -u^*- and the «i is the 
same except U[ nst ' and t^. One can show easily that the structures of the poles are the 
same between [J[ nst ' and U\, and hence the possible difference is a regular term. From 
the gauge invariance and the dimensional counting, this regular term is restricted in the 
following form: 

N c 

U^Nc, N f ) = U[ nst \N c , Nf) + C No 5 Nf ,2 Nc -2 + D Nc 5 Nf>2Nc £ aj. (29) 

i=i 



5 For Nf = 2N C case, one notices that the square inverses of the gauge coupling constant 1/g 2 in the 
microscopic theory and l/fflx appearing in the exact solution are different by a constant shift in general. 
This is not an inconsistency even for the scale invariant case, because finite renormalizations might exist 
in general. 
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In fact, C3, D s ^ for all the above proposed curves. 

These differences do not seem to lead to any inconsistencies because, for the Nf = 
2N C — 2, 2N C cases, the anomaly free symmetries certainly allow the above ambiguities in 



the construction of the curves |18|, [T|j, |2(|. There would also be some non-perturbative 
ambiguities in the definitions of the quantum operator TrA 2 in the microscopic theory^. 
For N c = 2 case, this type of resolution has been discussed in f22|| . We point out that 
these differences cause some qualitative differences in the non-perturbative renormaliza- 
tion procedures under reductions of the systems. Consider for example the Higgs breaking 
di = a'} — b (i = 1, • • • , N c — 1), Oj\r c = (N c — 1)6 (b — > 00) while K matters are kept mass- 
less; rrii = i\/2b (i = 1,---,K). Then the system with (N c ,Nf) reduces to that with 
(iV c — 1,K). When it is applied to the microscopic result (|16D, there appear divergences 
which cannot be absorbed into the matching condition (plf). But it can be easily shown 
that the divergences are regular terms and fllTf ) becomes consistent if the operator u is 
subtracted by divergent regular terms u™^ti k = u Nc iv — {regular terms) under the re- 
duction. On the other hand, we do not need such a subtraction for the modulus u under 
the similar reduction in the curves proposed in refs. [0, 20fl . Another example is the 
reduction in the mixed branch. If one sets a^ c = in the curve of Jl^] with (N c , Nf) 
for the massless case, the curve reduces to the one with (N c — l,Nf — 2). This can be 
justified physically by assuming that the mixed branch ajv c = with q l a = (g)^/-!^ > 
Qi = (q)^ 1 '^n c touches the Coulomb branch at (q) ~ ]19[. On the other hand, the 
microscopic one-instanton result is not consistent with naively setting ajv c = 0, and u 
needs again regular term shifts. A further analysis would be necessary in order to clarify 
the physical meaning of this difference. 



After completion of this work, we noticed the preprints p3| , in which the explicit eval- 
uations of the prepotentials of the exact solutions are discussed for iV = 2 supersymmetric 
gauge theories. 



6 One would have u mst - = c(q)u (c(0) = 1) and u lnst - — u + c'A 2 for massless Nf — 2N C and 
Nf = 2N C — 2 cases, respectively, which could explain the discrepancies in (|(|). For N c = 2 case, see 
refs. l2ll 
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